The SU(2) Gauge Theory of the Substantially Doped Mott-Insulator: 
the phase separated d+s superconductor. 
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We reformulate the mean-field theory associated with the SU (2) gauge theory of spin-charge sepa- 
ration in doped Mott-insulators of Wen and Lee so that it can deal with the regime of substantial 
carrier density. We find that it describes remarkably accurately phase separation tendencies. More- 
over, for elegant reasons the theory insists that the superconducting state should have a d + s 
orderparameter with a s-wave component growing in the overdoped regime. This appears to be 
consistent with both Raman- and tunneling measurements in cuprate superconductors. 
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The d-wave nature of the superconducting order pa- 
rameter in the cuprates is often taken as the leading evi- 
dence for this superconductivity to be caused by an non- 
phonon mechanism rooted in strong electron-electron in- 
teraction physics. However, a number of evidences ap- 
peared in independent experiments indicating that the 
order parameter breaks the symmetry further. The d- 
wave order parameter respects the rotational symmetry 
of the tetragonal (square) lattice, but by the admix- 
ing of a s-wave component this is spontaneously broken 
to an orthorombic symmetry. In orthorombic crystals 
like the YBCO superconductors such a s-wave admix- 
ture is automatic (and well established!]]]) but both in 
Raman- 2 and c-axis tunneling [3] experiments evidences 
have been found for a substantial (up to 20%) s-wave 
admixture to occur spontaneously in the BISCO super- 
conductors characterized by a square Cu — O lattice. 

It seems that in none of the various candidate theo- 
ries for the mechanism of superconductivity one finds a 
truly good reason for d + s symmetry. The 'slave' SU(2) 
gauge theory due to Wen, Lee and coworkers^ is among 
the most sophisticated unconventional high Tc theories. 
Motivated by the question of how to deal with charge in- 
homogeneity in this theoretical setting we found out how 
to compute in the regime of substantial doping, stay- 
ing strictly within the original formulation of the theory. 
On the one hand, we find that the ST/ (2) theory deals 
remarkably well with phase separation, but we also dis- 
covered that for elegant reasons rooted in the projective 
symmetry structure |5| behind the theory, the supercon- 
ductivity has to be d + s in the overdoped regime. Our 
main findings are summarized in the figures. As function 
of increasing J ft we find that a phase separation regime 
opens up at low doping (Fig. 1) saturating in a complete 
phase separation at J/t ~ 4, consistent with the out- 
comes of exact numerical diagonalization studies @, 0|- 
Also, the electronic incompressibility calculated for a re- 
alistic J/t = 0.3 is remarkably consistent with experi- 
ment (Fig. 2). This suggests that by adding physics 
beyond the t — J model (e.g., long range Coulomb in- 
teractions) the SU(2) theory can deal with the stripes 



USES EH- Above a doping x = 0.15- 0.20 a homoge- 
neous superconducting phase with d + s symmetry is pre- 
dicted, characterized by an s-wave component growing 
proportional to the doping and a d-wave component that 
is decreasing, to vanish completely at a doping x ~ 0.3 
(inset Fig.l). 

The above was overlooked in the earlier 

work[l[]3] 

for 

mere technical reasons: to simplify the calculations the 
focus was on the low hole density limit where the charge 
sector can be described as a dilute bose gas. However, 
at higher dopings one has to do justice to the hard-core 
nature of the bosonic charge sector and this is actually 
quite easy to implement on the mean-field level. Be- 
fore we go into the details of the calculations let us first 
present a simple argument making transparent why the 
SU(2) projective symmetry structure (i.e., the mean-field 
states span up the gauge volume [5]) implies that the su- 
perconducting order parameter has to be d+s. The idea 
to invoke an enlarged (SU(2)) gauge symmetry is rooted 
in the observation by Affleck and Marston[13( that the 
spin system at half-filling in a fermionic 'spinon' rep- 
resentation is characterized by both the usual 'stay at 
home' t/(l) gauge symmetry, and a local charge conju- 
gation symmetry meaning that one can describe the spin 
system equally well in terms of spinon particles and anti- 
particles. In this representation the spin system is in- 
variant under SU (2) gauge transformations tpi — > giipi 
where ipi = [fit, fh) (with p creating a spinon) and 
cji G SU(2). Away from half- filling the charge conju- 
gation is broken and SU(2) — * U[l). However, the 
St/ (2) symmetry is only weakly broken and the gauge 
transformations away from the 'equator of the isospin 
sphere' turn into soft physical modes [l2j]- To encode this 
mathematically, Wen and Lee introduced a doublet of 
bosonic holons hi — (6ii, &2i) transforming covariantly 
as hi — > Qihi. To appreciate the meaning of the boson 
SU(2) flavors it is informative to inspect how the phys- 
ical (gauge invariant) Hilbert space looks like[4[. The 
singly occupied states correspond with one spinon per 
site f\ |0), while the physical hole corresponds with the 
interesting combination -^m[o[ +^2/{/|) l^)- This reveals 
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the meaning of the boson doublet. The hole is a spin- 
singlet but the SU(2) parametrization remembers the 
charge conjugation symmetry of half-filling, insisting that 
singlets formed from missing spinons or of local s-wave 
pairs of spinons are gauge equivalent. The physical hole 
is a coherent superposition of these two possibilities and 
the holon flavor takes care that the state as a whole is a 
gauge singlet. 

On the mean-field level the spinon sector takes care 
of the symmetry of the superconducting order param- 
eter. Let us consider the completely empty state: 
the mean-field theory is of course unphysical because 
the gauge fluctuations will have destroyed the de- 
confining spin-charge separated state already at much 
lower hole densities. But it is an informative exer- 
cise. The exact wave function of the empty state is 

\vac) phys = IT \ b\x + b\ 2 fl l fl^j \vac) unphys . Spin-charge 
separation implies that the system looses its knowl- 
edge about the three particle correlation ~ b\ 2 f^f^. 
In addition, one has to satisfy the local constraint 
(iplT l ipi + h\T l h^ |phys) = 0,/ = 1,2,3, implying that 

the bosons have a hard core: at most one boson of either 
flavor can be present at a single site. The best choice for 
a holon-spinon product wave function becomes obvious: 
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I ( b li + b i2) ( l + fijii )\ \vac) unphys . 
In this mean field state there is exactly one boson present 
and it follows that this state decribes a holon Mott- 
insulator together with a local pair s-wave spinon super- 
conductor. Of course, this state is nonsensical ('nothing- 
ness' is not a superconductor) but it does remarkably well 
as a density functional theory 14| ground state. The en- 
ergy of the 'exact' empty state is of course vanishing and 
the same is true for the above mean field state: the hard- 
core of the bosons kills the hopping energy, and since the 
spin-operators have the form S l — fa<J l a pfp, these vanish 
as well. Lastly, it is readily checked that the three projec- 
tion operators tpj^tpi + h[T l hi also act like zero. In what 
follows we will merely use the slave mean-field theory as 
a density functional theory and the 'magic' revealed in 
the above gives us reasons to believe that it might work 
better than naively expected. 

Let us now turn to the full mean-field theory. The 
starting point is the usual t — J model and after the 
Hubbard-Stratonovich decoupling in the spinon sector in 

terms of the auxiliary fields Xij = {^f\ihj + f\%fij^ an d 

^ij = (f]iflj ~ fliflj) one obtains the standard result, 

H ™f = h i k i - fl 0* Q^tkrf + h W h ^ (!) 

E Y (J'-' V + 4Ui^o) + thtUijhj + h.c. 



<i]> 



and iUij E SU(2), up to a scalar factor. The mean- 
fields are invariant under the projective transformation 
Uij -> Uij = .'/.' .,.'/, _/ = gla l 0l T l gi. This implies 
that we may choose some reference state U^, and label 
all other mean-field states by gi. We choose as refer- 
ence mean-field state the d-wave superconducting state 
(dSC), Ul + - = -XT 3 + At 1 , U? i+ij = - X r 3 - At 1 . The 
mean- field solutions in the projective gauge volume [5[ 
are characterized by the isospin SU{2) sphere; we take 

g i — exp ^— i-^-T 2 ^j such that the fermionic latitude angle 

9 F = ±|tt decribes the staggered flux phases and 9 F = 
the superconductor, while the longitude (<£>) represents 
the 17(1) gauge[||. 

In Eq. (TTJ the holons are still fully dynamical. How- 
ever, we learned from the empty limit that their hard- 
core nature is crucial and it is quite simple to construct 
condensate wavefunctions for hard-core bosons. These are 
just like XY spins and the mean-field states have the gen- 
eralized coherent state form, taking care of the SU(2) 
doublet: 

\^o) holons = II f (oti + Piiuib^ + Uj&] 2 )J \vac) unphys (2) 

where itj = cos(^-),u, = sin(^), and \a\ 2 = 1 — |/3| 2 . 
Further, ft — > /3oie l ' Pi such that fth is real in terms 
isospin angles 6, <j>. For clarity's sake, we distinguish the 
bosonic angle from the fermionic isospin angle F . This 
parametrisation is convenient, since 6% = ^ makes the 
expectation values for b\ with vacancies indistinguish- 
able from 62 with a spinon pair, reproducing the particle- 
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level, the local constraint equation is replaced by its av- 
erage: (ipl^-ipi + h\r l hi^ = 0, and for the condensate 
Ansatz Eq. ([2]) this becomes 

= |/3 0l | 2 sin(^)cos(^) (3) 

= -i|/3 i| 2 sin(6' i )sin(^ i ) (4) 

= |/3 0j | 2 cos(^)- (5) 

Eq. (HJ) is satisfied by the f7 (1) gauge (pi = 0, such that 
= 0, and Eq. ((5]) is standard. The fermionic expecta- 
tion values are calculated with respect to the dSC gauge. 
However, in the original treatment by Wen and Lee[4|, as 
well as in the follow up work [3] the first constraint equa- 
tion Eq. © appeared to be automatically satisfied as a 
consequence of the treatment of the holons in terms of the 
dilute bose-gas limit . This is not at all true at substantial 
dopings, as can be immediately seen from Eq. ([3]): for 
any finite doping (|/?oi| 0); m the presence of a d-wave 
superconducting order parameter (6 = -|), one needs to 
have a finite s-wave component that is linearly propor- 
tional to the doping: < fLfL >— § ! The presence of 
d-wave superconductivity is of course not necessary at all 
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for the s-wave superconductor, as is immediately obvious 
from the empty limit. 

Let us now explicitly derive the free energy functional 
from the Hamiltonian fTJ). We have checked that at least 
for the t — J model homogeneous states are preferred and 
here we specialize to 0, = 9, \/3 ni \ = [3 and a l 0i — a l , while 
Xij and Aij are chosen to describe a homogeneous flux 
phase and/or d-wave superconductors. The ground state 
energy per site &mf becomes 
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E k + —J(\ X \ 2 + \A\ 2 ) 



(6) 



- (/i + ao sin0 + aocos60|/3| 2 



Here E(k) = \J (xk — a^) 2 + (A fe — aj) 2 is the spinon 
dispersion relation, with Xk = ^-x( cos ^x -t-cosfc,,) and 
Afe = ^-A(cosk x — cosky). Notice that we consider the 
mean-field energy density functional in the grand canon- 
ical ensemble, by fixing p instead of the density p, antic- 
ipating phase separation. Using the notation p = |/3| 2 , 
the four saddle point equations are 
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, . s 1 / u + aj sin 9 + cos ( 
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where the third equation is the homogeneous form of Eq. 
([3|). We find solutions by minimizing the energy ([6]) nu- 
merically with the simulated annealing method [Hsf . We 
also calculated the saddle point solutions for a staggered 
ajj, to incorporate the possibility of the staggered flux 
state 13]. Due to Eq. ([3]), a staggered leads to a 
staggered s-wave component. 

Some main features survive in this revised version of 
the mean-field theory: in the presence of holes the su- 
perconducting state is always preferred over the flux 
phased 

, and at least for the t — J model it appears 
that uniform states have always a better energy than in- 
homogeneous states 12j . A first novelty is that we now 
find that macroscopic phase separation occurs generically 
at low hole density: as function of chemical potential, the 
system stays initially at half-filling and pending the ratio 
of J/t at some finite p a level crossing takes place to a 
state with a finite doping level (Fig.l ). As function of 
increasing J/t the width of this phase separation regime 
is increasing (see Fig. 1) and we find that for J/t oi 4 the 
phase separation is complete, consistent with exact diag- 
onalization studies on the t-J model indicating a com- 
plete phase separation for J/t > 3.5 [H, 0]. This is by 
itself a rather serious test of this mean-field theory in its 
capacity of a density functional theory: to reproduce the 
critical J/t for complete phase separation, both the ener- 
gies of the Mott-insulator (x = 0) and of the states having 
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FIG. 1: The zero temperature phase diagram for the SU(2)- 
mean field theory for the t — J model, as function of doping 
x versus the ratio of exchange (J) and hopping (t) energies. 
There is a regime to the left of the blue line where the system 
macroscopically phase separates in the Mott-insulator and a 
homogeneous d+s wave superconductor. This d+s supercon- 
ductor only exists in a rather small doping regime adjacent to 
the phase separation line, and a transition follows to a pure 
s-wave superconductor at higher doping (green line, see the 
inset for the doping dependence of the gaps Ad, s for a 'phys- 
ical' J/t — 0.3). The false colors code for the overall gap 
magnitude, and these together with the d and s gaps in the 
inset are computed in the canonical ensemble in the phase 
separation regime, to indicate how they would evolve when 
phase separation would not take place. 



a hole density i ^ 0,1 have to be quite accurate, given 
that the mean-field energy of the empty state (x = 1) is 
exact. How well does this compare to experiment? In 
Fig. 2 we compare the calculated doping dependence of 
the chemical potential for a 'physical' J/t = 0.3 with 
the photoemission results by Fujimori and coworkers [18j. 
The doping independence of p is indicative of phase sep- 
aration and the theory predicts a wider phase separation 
regime (x cr it ~ 0.2) than experiment {x cr u ~ 1/8): this 
discrepancy has likely more to do with the t — J model 
itself than with the mean-field theory. However, we find 
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that the electronic incompressibility 1/k — — 
of the homogeneous superconductor at higher dopings is 
remarkably well reproduced by the theory (Fig. 2). 

Let us now address the nature of the order parame- 
ter of the homogeneous superconductors found at higher 
dopings (inset Fig. 1). We already established that the 
s-wave component < fLfh >= f , as a direct conse- 
quence of the constraint Eq. ([3]). To further emphasize 
this matter, we compare in the inset of Fig. 2 the energy 
of a state where we have fixed the Lagrange multiplier 
= such that the s-wave component vanishes, with 
the best d + s state, finding that the former is indeed 
a false vacuum. The surprise is that according to the 
SU(2) mean-field theory a phase transition is present at 
a higher doping where the d-wave component vanishes 
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FIG. 2: The chemical potential shift as a function of dop- 
ing computed for J/t = 0.1,0.3. Using t = 0.44 eV [H the 
absolute values of Afj, can be compared with the experimen- 
tal (photoemission) results by Fujimori et al taking the 
'physical' J/t = 0.3 the theory overestimates the critical dop- 
ing where phase separation stops, but the theory is accurately 
reproducing the electronic incompressibility (the slope of the 
lines) of the homogeneous superconductor. The inset shows 
the relative mean-field energy cost of both the uniform pure d 
wave superconductor (red) and the superconductor with stag- 
gered s-wave component (blue), relative to the homogeneous 
d + s superconducting ground state as a function of chemical 
potential. The spikes in the relative energy differences are due 
to the fact that for low chemical potential, the pure d-wave 
and the staggered d + s-wave superconductors do not admit 
a non-zero holon density. 



altogether. As can be seen from Fig. 1, a d-wave compo- 
nent is only found in a relatively narrow doping region 
adjacent to the phase separation line, regardless the J/t 
ratio. Notice that we show in Fig. 1 the development of 
the d and s wave components and the overall gap magni- 
tude as computed in the canonical ensemble: because we 
impose an average density these are actually unphysical 
in the phase-separation regime. 

Summarizing, by just improving the technical formu- 
lation of the SU(2) gauge theory on the mean-field level 
(holons as hard-core bosons) we have discovered that 
this theory exerts some remarkable, unexpected powers. 
Viewed as a density functional theory [14j, it is quanti- 
tatively accurate, even in a regime where its literal in- 
terpretation as a dynamical theory becomes nonsensical 
(the confinement problem at large doping). It reproduces 
the macroscopic phase separation tendency of the t — J 
model accurately and it is well understood that this is 
an unphysical pathology of the this model itself. It just 
signals the presence of stripe- like instabilities [1,0, [HI II | 
and we leave it for future work to find out what the gauge 
theory has to tell about stripes. 

The punchline is that the theory produces a strong- 
and counterintuitive prediction that is well in reach of 
experimental confirmation: for deep and elegant rea- 
sons specific for the workings of the underlying projective 
gauge principle the superconducting order parameter 
in the overdoped regime has to be d + s, with an s-wave 
component increasing linearly with doping and a decreas- 



ing d-wave component, that is governed by the proximity 
of a transition to a pure s-wave state at a doping level 
which is likely beyond the level that can be experimen- 
tally realized. This is surely not falsified by the present 
state of experimental affairs. We predict a substantial 
s-wave admixture at optimal doping of order 10 — 20%, 
which is consistent with c-axis tunneling experiments 0] . 
We predict that the gap ratio r = A^/A^ grows with dop- 
ing, in accord with the Raman measurements 0]. There 
are quite a number of other, in principle straightforward 
ways to measure the evolution of the d + s order param- 
eter and we hope that the opportunity to (dis)prove the 
Wen-Lee gauge theory will be sufficient motivation for 
such an experimental effort. 
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